Depletion potential in colloidal mixtures of hard spheres and platelets 
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The depletion potential between two hard spheres in a solvent of thin hard disclike platelets is 
investigated by using either the Derjaguin approximation or density functional theory. Particular 
attention is paid to the density dependence of the depletion potential. A second-order virial ap- 
proximation is applied, which yields nearly exact results for the bulk properties of the hard-platelet 
fluid at densities two times smaller than the density of the isotropic fluid at isotropic-nematic phase 
coexistence. As the platelet density increases, the attractive primary minimum of the depletion po- 
tential deepens and an additional small repulsive barrier at larger sphere separations develops. Upon 
decreasing the ratio of the radius of the spheres and the platelets, the primary minimum diminishes 
and the position of the small repulsive barrier shifts to smaller values of the sphere separation. 

PACS numbers: 61.20.-p, 61.20.Gy, 82.70.Dd 



I. INTRODUCTION 

Depletion interactions between big colloidal particles 
induced by smaller particles, which can be either the sol- 
vent particles or a colloidal component of its own, are of 
significant current research interest because of the impor- 
tance of these effective interactions in various colloidal 
processes. For example, flocculation of colloids can be 
driven by the addition of non-adsorbing polymers via the 
depletion mechanism pj . Whereas experimental and the- 
oretical studies have focussed on binary hard sphere flu- 
ids as well as on colloidal mixtures of hard spheres and 
hard rods or polymers, less attention has been paid to 
hard platelets acting as depletants, despite the great im- 
portance of colloidal platelets such as blood-platelets and 
clay minerals in both biomedicine and geophysics. Very 
recently a colloidal mixture of silica spheres and silica 
coated gibbsite platelets has been stabilized for the first 
time 0, and the depletion potential due to the presence 
of thin hard platelets has been derived theoretically for 
non-interactingplatelets corresponding to the limit of in- 
finite dilution jJQ- It has been found that the Derjaguin 
approximation for the depletion potential yields accurate 
results for non-interacting platelets provided the ratio of 
the radius of the spheres and the platelets is large [2|. 

In this paper we focus on the depletion interaction in- 
duced by thin hard platelets, taking into account the 
steric interactions between the platelets in terms of a 
second-order virial approximation. On the basis of our 
recent theoretical studies on fluids of thin hard platelets 
near hard walls 0, @ , we expect that excluded volume 
interactions between the platelets influence the deple- 
tion interaction already at rather low platelet densities 
due to their cumbrous shape as compared with spherical 
or rodlike depletants. Taking excluded volume interac- 
tions into account is particularly interesting because cor- 
relation effects may cause repulsive features of depletion 
forces which are important in the context of colloidal sta- 
bility In the present paper we use density functional 



theory (Sec. II) to study the depletion potential between 
two hard spheres induced by thin hard platelets (Sec. 
III). Particularly, we compare the results with the ones 
obtained for non-interacting platelets. 

II. DENSITY FUNCTIONAL THEORY 

We consider an inhomogeneous fluid consisting of thin 
platelets of radius R p in a container of volume V. The 
platelets are taken to be hard discs without additional 
attractive or repulsive interactions. The number density 
of the centers of mass of the platelets at a point r with 
an orientation uj — (9, </>) of the normal of the platelets 
is denoted by p(r,uj). The equilibrium density profile 
of the inhomogeneous liquid under the influence of an 
external potential V{v, uj) minimizes the grand potential 
functional 

0[p(r,w)] = J dr 3 dujp{r,uj) [k B T (ln[47rA 3 p(r, uj)] - l) 
- fx + V(r,cj)}+F ex [p{r,cj)}, 

(1) 

where A is the thermal de Broglie wavelength and p 
is the chemical potential. The free energy functional 
F ex \p(v, uj)] in excess of the ideal gas contribution has 
not been taken into account inprevious studies on the 
depletion force due to platelets 0, Lj ■ We express the ex- 
cess free energy functional as an integral over all possible 
configurations of two platelets 

F ex [p(r,oj)] = -~y~ J drldu)xdrldw2p{Tx,wi) 

x/pp(ri2,wi,w 2 )p(r2,W2) , (2) 

where 1*12 = ri— V2 and / pp (ri 2 , uj\, LJ2) is the Mayer func- 
tion of the interaction potential between two platelets. 
The Mayer function equals —I if the platelets overlap 
and is zero otherwise. Explicit expressions of the Mayer 
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function for thin platelets are documented in Refs. Q 
and pfl. 

For the homogeneous and isotropic bulk fluid the grand 
potential functional [Eq. |JIJ] reduces to 

^ = p b [k B T(H\ 3 Pb ] - 1) - m] + Y R3 P Pbk B T , (3) 

where p b — V^ 1 J dr 3 J du>p(r,u)) is the total particle 
number density. The equation of state derived from the 
grand potential [Eq. (JSJ] takes the following form: 



(4) 



with p\ — pbRf, and p\ = pbRp/iksT). The same equa- 
tion without the second term in parenthesis holds for 
the ideal gas limit (i.e., non-interacting platelets). With 
increasing particle number density the ideal gas equa- 
tion of state on one side and the second-term virial series 
[Eq. as well as computer simulation data 0, H, Q on 
the other side deviate. Thus for p\ = 0.1 the osmotic 
pressure pi = 0.15, as calculated from Eq. J3J|, agrees 
exactly with simulation data, while the ideal gas equa- 
tion of state underestimates the osmotic pressure by a 
factor of 1.5. The comparison of the calculated equation 
of state with computer simulation data exhibits that the 
two-term series in Eq. Q is a good approximation for 
pi < 0.2, whereas the ideal gas model may be used for 
very low particle number densities p* h < 0.04. For a dis- 
cussion of higher order virial terms for fluids consis ting 
of hard platelets we refer to Refs. and [13, IH Ell- 
in the present study we restrict our attention to parti- 
cle number densities p\ < 0.2 for which the second-order 
virial approximation is appropriate and the platelet fluid 
is in the isotropic phase. For comparison, the isotropic- 
nematic phase transition is first order with coexistence 
densities pbiR 3 — 0.46 and PbNRp = 0.5 according to a 
computer simulation |9(. 



III. THE PLATELET-INDUCED DEPLETION 
POTENTIAL BETWEEN TWO SPHERES 

The results of the preceding section show that inter- 
molecular interactions between platelets increase the os- 
motic pressure of the bulk fluid already at low particle 
densities. Now we study the influence of intermolecular 
interactions on the depletion potential between two hard 
spheres of radius R s immersed in a fluid of hard platelets 
of radius R p . The depletion potential W{h) is the free 
energy difference between the configurations of two big 
spheres at fixed distance h immersed in the solvent and 
at macroscopic separation h = oo (see Fig.^l. 



A. The Derjaguin approximation 

The depletion potential W(h) between two hard 
spheres at close distance due to the presence of small 





FIG. 1: The system under consideration consists of two hard 
spheres of radius R s immersed in a solvent of hard platelets of 
radius R p . The separation between the surfaces of the spheres 
is denoted by h. Only the projection of the spheres on the 
plane of the figure is shown. 



platelets (R p C fl s ) can be calculated from the finite 
size contribution of the grand potential function uj(h') 
of the platelet fluid confined between two parallel hard 
walls at distance hi using the Derjaguin approximation 




W, 



Derj \ 



(h) =ttR s / dh'uj(h'). 



(5) 



where h is the separation between the surfaces of the 
spheres. (For the subtle issue of the range of validity of 
the Derjaguin approximation see Refs. |Q and [Hf-i 

We first consider a hard-platelet fluid confined by two 
parallel hard walls at z = and z = h, and calculate the 
surface and finite size contributions to the grand poten- 
tial defined via 



n[p(z,e,(f>)} = Viub + 2A 7 + Aw{h) , 



(6) 



where A is the area of a single surface, LOb is the bulk 
grand canonical potential density, and V is defined as 
the volume of the container with its surface given by the 
position of the rim of the particles at closest approach 
so that V — Ah. 7 is the wall-liquid surface tension in 
the absence of the second wall and ui(h) is the finite size 
contribution. Figure[2]displays the calculated surface and 
finite size contributions to the grand potential together 
with the results for non-interacting platelets 0, El 



[27 + "WU„ t 




0<h<2R v 



h > 2R ri 



(7) 

As is apparent from Fig. [3 the steric interaction between 
the platelets increases the surface contributions with in- 
creasing density. Within our numerical precision, we 
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FIG. 2: The surface and finite size contribution 27 + ui(h) 
to the grand potential as obtained from Eqs. Q, @, and 
© [solid lines] of a fluid consisting of thin hard platelets of 
radius R p confined in a slit of width h and in contact with 
an isotropic bulk reservoir at density pt- The dashed lines 
represent the corresponding results for an ideal gas of platelets 
[see Eq. Q]. The width of the slit increases from bottom to 
top: h/R p = 0.5, 1.5. 



found that an accurate evaluation of the wall-liquid sur- 
face tension 7 could be achieved for a fixed wall sepa- 
ration h = AR p at all considered densities PbRp < 0.2. 
For much higher densities, a larger value of h might be 
required because of the wetting of the wall-isotropic liq- 
uid interface by a nematic film of diverging thickness 
On the other hand the wall-liquid surface tension 
lideail (ksT) = pbRpit /4 for non- interacting platelets fol- 
lows from Eq. J7J) for a wall separation h = 2R p . For a 
detailed discussion of the surface tension and the excess 
coverage as well as the density and orientational order pa- 
rameter profiles of fluids consisting of thin hard platelets 
near a single hard wall we refer to Ref. [4]. The results 
for the finite-size contribution u>(h) are shown in Fig. [31 
(a). As function of h the finite size contribution corre- 
sponds to the solvation free energy for the immersed two 
plates acting as the confining walls for the fluid and, by 
construction, uj(0) = —2-f and w(oo) = 0. Upon increas- 
ing the platelet density, the attractive minimum of uj(h) 
at h = deepens and a maximum at larger values of 
h develops. The corresponding solvation force per unit 
area f(h) = —dui(h)/dh is attractive for small slit widths 
h as is shown in Fig. [31(b). Upon increasing the platelet 
density, the cusp of the solvation force at h = 2R p sharp- 
ens. For comparison we note that the maximum at the 
cusp is more pronounced for the confined platelet fluid 
than for a corresponding rod fluid [6| due to the rela- 
tively larger steric interactions between platelets as com- 
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FIG. 3: (a) The finite size contribution u)(h) to the grand 
potential as obtained from Eqs. Q, ||2J, and (JSJ of a fluid 
consisting of thin hard platelets of radius R p confined in a slit 
of width h for three values of the density of the isotropic bulk 
reservoir: pi = pbRp = 0.05 (dotted line); pi =0.1 (dashed 
line); pi = 0.2 (solid line), (b) The solvation force per unit 
area f(h) = —du{h)/dh of the same fluid [with the same line 
code as in (a)] as a function of h. 



pared with those between rods. Moreover f(h) is a con- 
vex function for slit widths smaller than two times the 
radius of the platelets, while the solvation force in a sol- 
vent of hard rods is a concave function for slit widths 
smaller than the length of the rods [f| . Figure 01 displays 
the depletion potential together with the results for non- 
interacting platelets which can be calculated analytically 
from Eqs. © and 10: 
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FIG. 4: Depletion potential WDerj{h) between two hard 
spheres of radius R 3 = 3.5R P due to the presence of thin 
hard platelets of radius Rp as obtained from the Derjaguin 
approximation [Eq. @] for various bulk densities pb of the 
platelets: pi = pbR% = 0.02 (top curves); pi = 0.08 (middle 
curves); pi = 0.2 (bottom curves). The solid lines (with dots 
marking the maxima) represent the calculations for interact- 
ing platelets and the dashed lines denote the results for non- 
interacting platelets. For reasons of clarity, the lower four 
curves are shifted down by -0.3 fcsT and -0.9 fcsT, respec- 
tively. The dashed lines are zero for h/R p > 2. Although the 
interactions have only minor influence on the depth of the pri- 
mary minimum, the depletion potential becomes significantly 
less attractive for increasing densities p p . This weakens the 
platelet induced flocculation of a solution of big spheres. 



energy per sphere is 

f{ps, Pb ) = f (hs) {Ps) 

OO 

+2tt Ps [ drr 2 W(r-2R s )g^(r,p s ), (9) 

where g^ hs ^> (r, p s ) is the radial distribution function of the 
pure sphere fluid and f^ hs \ps) is the Helmholtz free en- 
ergy per sphere for the same homogeneous fluid of density 
p s . In the limit R, 



3> R p the radial distribution function 



(r, p s ) is almost constant over the range of integra- 



g (hs) 

tion where W(r — 2R S ) ^ and we can approximate it 
by its constant contact value g^ hs ^(2R B , p s ) [lg. Using 
Eq. JSJ resulting the integral in Eq. © can be evaluated 
analytically for non-interacting platelets: 
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The integrated strength of the depletion potential Apf^j 
is negative, reflecting the fact the depletion potential is 
always attractive (Fig. ■ A numerical calculation of the 
corresponding quantity A^erj for interacting platelets ex- 
hibits that the steric interacting between the platelets 
weakens the integrated strength of the depletion poten- 
tial by 22 % at a bulk density PbRp = 0.2. Hence the 
influence of steric interactions between platelets might 
be quite visible for phase equilibria. For example, the 
thermodynamic onset of flocculation of colloidal spheres 
induced by the depletion effect will be reduced due to 
platelet interactions. 



ing platelets exhibits a small barrier at larger sphere 
separations h in addition to the primary minimum at 
h = 0. With increasing platelet density the depletion po- 
tential deepens and the position of the maximum shifts to 
smaller values of h. The small repulsive barrier will have 
minor effects on kinetic stabilization, although the repul- 
sive features might still be measurable. For example, our 
numerical calculations exhibit a maximum barrier height 
of 0.25 ksT relative to zero at a density of PbRp = 0.2 
in a system of size ratio R s /R p = 3.5 corresponding to 
the aforementioned mixture of silica spheres and gibb- 
site platelets 0- In order to examine the influence of 
steric interactions between platelets on thermodynamic 
properties of sphere-platelet mixtures, we treat the de- 
pletion potential as a perturbation to the hard-sphere 
potential. The first-order approximation in this thermo- 
dynamic perturbation approach for the Helmholtz free 



B. Density functional approach 

A general approach for calculating the depletion po- 
tential is based on a density functional theory (DFT) for 
a mixture of hard spheres and the particles acting as de- 
pletants [H [3. This approach avoids the Derjaguin 
approximation. Translating the particle insertion idea 
developed in Refs. and to the present system 
leads to the following expression for the depletion po- 
tential between two spheres due to the presence of hard 
platelets: 

W(r) 1 f 

-y^7=^ / dr x dui \p(ri,u)) - p(oo,w)]f ap (T-T U u) 

(11) 

where / sp (r — Ti,u>) is the Mayer function of the inter- 
action potential between a sphere and a platelet. The 
Mayer function equals —1 if the particles intersect or 
touch each other and is zero otherwise. p(ri,u>) is the 
density profile of platelets in the external potential of 



a single fixed hard sphere located at the origin of the 
coordinate system and p(oo, ui) is the corresponding den- 
sity profile of the bulk fluid. We emphasize that the 
density profile entering Eq. (|ll|l depends only on equi- 
librium properties of the depletant fluid in the absence 
of the second hard sphere to be inserted at position r. 
This observation simplifies the calculation of W(r) con- 
siderably, because the symmetry of the density profile is 
determined solely by the symmetry of the external poten- 
tial of a single sphere fixed at the origin of the coordinate 
system. 

Apart from possible surface freezing at high densities, 
non-uniformities of the density depend only on the ra- 
dial distance r = |r|, so that p(r,iv) — p(r,u>). Hence, 
calculating density profiles before insertion of the second 
hard sphere is much easier than after insertion, when the 
presence of the second sphere leads to a more complex 
spatial variation of the densities. A detailed discussion 
of Eq. 1)1 1H and its application to the analogous case of 
non-interacting hard rods acting as depletants is given in 
Ref. 0. 

For an ideal gas of platelets in contact with a fixed hard 
sphere the density profile reduces to p(ri , u>) — p(oo, u>) = 
f sp (r 1: uj) so that the integral in Eq. ((TTJl has a purely 
geometrical meaning and measures the excluded volume 
of a platelet confined between two hard spheres located 
at the origin of the coordinate system and at position r, 
respectively. 

In order to take intermolecular interactions between 
the platelets into account we first calculate numerically 
the density profile p(r, u>) of platelets in an external po- 
tential of one fixed hard sphere of radius R s . There- 
after the integral in Eq. Ijllfl is evaluated by inserting 
this density profile. To our knowledge, this technique 
has not been used before for interacting non-spherical 
colloids acting as depletants. 

The orientational averaged density profile 

p(r) = J <kjp{r,u) (12) 

of the platelet fluid in contact with one fixed hard sphere 
is shown in Fig.JSJfor various radii R s of the sphere. Upon 
increasing r > R s from the surface of the sphere the av- 
eraged number density increases and exhibits a cusp at 
r = R s + R p where platelets with their normal perpen- 
dicular to the radial direction touch the surface of the 
sphere with the rim. The maximum at the cusp is about 
25% above the bulk value pbRp = 0.085 for a size ratio 
R s / Rp — 5 and is less pronounced for smaller size ratios. 
The averaged density close to the surface of the sphere 
is larger for a small sphere than for a big one. Figure 
[fj] displays the calculated depletion potential for two size 
ratios R s /R p as a function of the separation between the 
surfaces of the spheres h = r — 2R S . With decreasing size 
ratio the range and the depth of the primary minimum 
shrinks and the position of the small repulsive barrier 
observed at higher densities shifts to smaller values of h. 
Moreover, the height of the repulsive barrier decreases 
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FIG. 5: Orientationally averaged density profile p(r) as ob- 
tained from Eq. W'll for hard platelets of radius R p in contact 
with a single hard sphere of radius R 3 located at r = 0. The 
radius of the sphere increases from left to right: R 3 /R p = 
0.5, 1, 3, 5. The arrows mark the location of the surface of the 
sphere at r — R s . Since the platelets are arbitrarily thin their 
density is nonzero for r > Rs- All curves exhibit a cusp at 
r — R s + R p followed by a decay towards the bulk density 
pbR? = 0.085, which is essentially reached at r = R s + 2R P . 
When the center of a platelet is located less than R p from the 
sphere surface, there are fewer possible orientations available 
to the platelet. 



upon decreasing the size ratio. These results are due to 
the fact that the number of platelets contributing to the 
depletion potential decreases as the ratio of the radius of 
the spheres and the platelets becomes smaller at a fixed 
bulk density (see Fig - 

In agreement with a recent theoretical study based 
on an evaluation of the excluded volume of a single 
platelet confined between two hard spheres, we find that 
the Derjaguin approximation for the depletion potential 
in the presence of non-interacting platelets yields accu- 
rate results for large size ratios R s /R p > 1. However, 
there are substantial deviations at higher densities as can 
be seen from Fig.[fj](a). The absolute value \W(h = 0)| 
of the DFT solution at contact is smaller than the one 
obtained from the Derjaguin approximation, and the re- 
pulsive barrier is less pronounced. Hence, for a size ra- 
tio Rg/Rp = 3 and a platelet density pbRp = 0.2 the 
Derjaguin approximation overestimates the depth of the 
depletion potential at contact by 37 % and the height of 
the repulsive barrier by 23 %. Increasing the size ratio 
R s / Rp > 3 does not lead to a significant improvement of 
the Derjaguin approximation as compared to the results 
shown in Fig. |H|(a). For R s /R p < 1 and higher densi- 
ties pb the DFT results deviate strongly from the predic- 
tions for non-interacting platelets [see p\ = pbRp = 0.2 
in Fig. ©(b)] . 
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FIG. 6: Depletion potential W(h) between two hard spheres 
of radius R a = 3R P in (a) and R s — 0.5R P in (b) due 
to the presence of thin hard platelets of radius R p as ob- 
tained from Eq. (I1H . The bulk densities of the platelets are: 
pi = PbRp = 0.02 (top curves); pi = 0.08 (middle curves); 
pi — 0.2 (bottom curves). Here h — r — 2R S is the sepa- 
ration between the surfaces of the spheres. The solid and 
dashed lines represent the calculations for interacting (DFT) 
and non-interacting platelets (ideal gas), respectively. In ad- 
dition the depletion potential as obtained from the Derjaguin 
approximation [Eq. JKJ] for interacting platelets is displayed 
by dotted lines in (a). For the smallest bulk density the solid 
and dotted line nearly coincide in (a). Since in (b) the spheres 
are only half as big as the platelets, in this case the Derjaguin 
approximation is unsuitable and therefore not shown. Only 
for sufficiently large densities a maximum denoted by a dot 
occurs. For reasons of clarity, the lower sets of curves are 
shifted down by -0.3 ksT and -0.9 ksT, respectively in (a) 
and by -0.05 ksT and -0.15 fcsT, respectively in (b). 



IV. SUMMARY 

We have applied a density functional theory to fluids 
consisting of thin hard platelets confined between two 
hard spheres (Fig. . Within the framework of a second- 
order virial approximation of the excess free energy func- 
tional, the depletion potential between the two spheres 
due to the presence of the platelets is determined numer- 
ically and compared with the corresponding results for 
non-interacting platelets. The main conclusions which 
emerge from our study are as follows. 

(1) Figure [21 demonstrate that steric interactions be- 
tween thin platelets of radius R p confined between two 
parallel hard walls increase the sum of the surface and fi- 
nite size contribution to the grand potential significantly 
already at rather low platelet densities PbRp ^ 0.025. 

(2) As function of the slit width h the finite size con- 
tribution to the grand potential of a slap of platelets ex- 
hibits a minimum at h = [Fig. El (a)]. A maximum at 
larger values of h is found for higher platelet densities. 
The corresponding solvation force is attractive for small 
slit widths and exhibits a cusp at h = 2R p [Fig. (b)] . 

(3) The depletion potential between two spheres as cal- 
culated from the Derjaguin approximation exhibits an 
attractive primary minimum at contact which deepens 
upon increasing the platelet density. Moreover, a small 
repulsive barrier at larger sphere separations develops 
with increasing density (Fig.^J. We find that the deple- 
tion barrier relative to zero is typically less than the ther- 
mal energy fceT, and therefore unlikely to significantly 
alter the kinetics of aggregation of the hard spheres at 
platelet densities smaller than two times the density of 
the isotropic phase at bulk isotropic-nematic coexistence. 
Nonetheless, with increasing platelet density the inte- 
grated strength of the effective interaction between the 
spheres becomes significantly weaker and thus reduces 
the thermodynamic onset of flocculation. 

(4) The orientational averaged density profile of a 
platelet fluid in contact with a single fixed hard sphere 
decreases towards the surface of the sphere because the 
range of accessible orientations is reduced when the par- 
ticle approaches the sphere. It exhibits a cusp at the 
position where the platelets lose contact with the sur- 
face of the sphere (Fig. EJ) . The maximum at the cusp 
decreases and the averaged density close to the surface 
increases as the ratio of the radius of the sphere and the 
platelets becomes smaller. 

(5) With decreasing ratio of the radius of the spheres 
and the platelets, the primary minimum at contact and 
the small repulsive barrier of the depletion potential di- 
minish and the position of repulsive barrier shifts to 
smaller values of the separation of the surfaces of the 
spheres (Fig. ISJ. From our numerical results based on 
a density functional theory for a mixture of spheres and 
platelets we found that the Derjaguin approximation is 
valid for large size ratio and very small platelet den- 
sity, but there are substantial deviations from the density 
functional results at higher densities even for large size ra- 
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tios [Fig. 0(a)] . These deviations increase with increasing densities the ideal gas approximation for the platelets 
platelet density. For small size ratios and high platelet becomes unsuitable [Fig. 0(b), p£ = 0.2]. 
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